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1. Introduction 

m 

The effect of spin-orbit coupling is well known from elementary quantum mechanics and 
Q\ '. atomic physics: it arises from the interaction between the spin of an electron in an atom 

and the magnetic field created by the orbiting nucleus in the rest frame of the electron. 
OS ! This magnetic field is related to the electric field created by the nucleus in its rest frame. 
If this field is a spherical electrostatic field, the interaction hamiltonian is given by 

Here, as usual, S and L are the spin and orbital momentum operators for the electron, 
m and e stand for the electron charge and mass, c is the speed of light in the vacuum 
and V(r) is the electrostatic potential of the atomic nucleus. For one-electron atom, 
the formula (|I|) is exact, otherwise V(r) can be thought as an approximation to an 
average radial potential experienced by the electron. Equation (|I|) is obtained in the 
non-relativistic limit (electron velocity is small compared to c — see, for instance, Q ) 
and so it is used in the non-relativistic description of an electron, i.e., by adding it to 
the Hamiltonian in the Schrodinger equation. 

In this paper we propose to examine a similar coupling that arises due to the 
relativistic treatment of the electron (i.e. using the Dirac equation) even in the absence 
of an external field. This is a consequence of the linearity of the Dirac equation in 
the space derivatives (and thus in the linear momentum operator p) and from the 
related fact that one needs a 4-component spinor to describe the electron. We will 
make the relativistic nature of this coupling apparent by solving the Dirac equation 



On the relativistic L — S coupling 



2 



in an infinite spherical potential well. Although the particle motion inside the well is 
free, the relativistic L — S coupling exists and vanishes only in the non-relativistic limit, 
which we are able to approach continuously by varying the well radius. In this limit the 
two-component spinor description is valid. 

A comparison between relativistic and non-relativistic solutions was already studied 
in |2j for a one-dimensional infinite square well potential. In the present paper we use 
the same procedures as in || to provide a bridge between known relativistic and non- 
relativistic solutions in the 3-dimensional spherical case, with special emphasis on the 
L — S coupling. Berry and Mondragon || have also applied similar methods in the 
framework of the Dirac equation in two spatial dimensions. 

In section 2 we pedagogically review the solutions of the free Dirac equation 
with spherical symmetry, in a slightly different fashion from the usual treatments, 
emphasizing the role of the L — S coupling term and its consequences for the set of 
quantum numbers of the solution. In section 3 we solve the Dirac equation for a spherical 
potential well and compare it to the non-relativistic solution of the corresponding 
Schrodinger equation for several well radii. Technical details, included for completeness, 
are mostly left to Appendices. 



2. Solutions of the free Dirac equation with spherical symmetry 

The free Dirac equation for a spin-| particle with mass m is a matrix equation for 
4-component spinors \& given by 

ih— = a ■ pc^ + f3mc 2 ^ (2) 

where p = — i ftV is the linear momentum operator, and a and (3, in the usual 
representation, are the 4x4 matrices 

"-(it) >-(;-,)■ 



Here / is the 2x2 unit matrix and a denotes the three Pauli matrices <7j i = 1,2,3 
obeying the relations 

a i°j = $ij + i £ ijkO-k i, j = 1, 2, 3 (4) 

where e^k is the anti-symmetric Levi-Civita tensor (£123 = 1) and summation over 
repeated indexes is implied. 

Using (£|) we can obtain the following general property of the a matrices 

a- A a- B = A- B + iAx B-t , (5) 

where A and B are two arbitrary vectors whose components commute with the matrices 
CKj and 
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is the 4-dimensional analog of the Pauli matrices. Using (|5]) and a ■ f a ■ f = I (here, of 
course, I stands for the 4x4 unit matrix), one can write 

a ■ p — a ■ r a ■ r a ■ p 

= a- r(f-p + irxp-'E) 

= a ■ f {f ■ p H — L ■ E) , (6) 

where r — \r\, f — f/r and L = f x p is the orbital angular momentum operator. 
Inserting (||) into the Dirac equation (H) we get 

i 

ih— = a ■ f (f -p + -L ■ E) c\I> + (3mc 2 ^ . (7) 
ot r 

Since the spin angular momentum operator in the Dirac formalism is S = | E the last 
expression contains a term involving the dot product L ■ S as in the spin-orbit term 
(H). This term is responsible for the L — S coupling^] for a relativistic spin-| particle 
even in the absence of an external potential. Clearly this fact is connected to the 
spinor structure of the wave function for spin-| particles and to the linearity of the 
Dirac equation, which leads to the appearance of the term a ■ p. In the Klein-Gordon 
equation, which is quadratic in the space derivatives, there is no such coupling. 

We can gain further insight into the origin of this effect if we write the spinor \1/ in 
(0) as a set of two-component spinors x' an d <p'- 




From the block off-diagonal form of the «j matrices one sees that the term a ■ p mixes 
the spinors x' an d V 9 '- m the literature (see, for instance, Bjorken and Drell 0]) this 
kind of term is referred to as "odd" as opposed to the terms like f3mc 2 which are called 
"even" , since they don't mix upper and lower two-component spinors. This distinction 
is important when one goes to the non-relativistic limit of the Dirac equation using 
the Foldy-Wouthuysen transformation [|], [5[ which aims at eliminating the odd terms 
through a unitary transformation and so decouple the upper and lower spinors. In this 
way one can regain the (Pauli) non-relativistic description of a spin-^ particle. The fact 
that the term responsible for the L — S coupling is odd indicates that this relativistic 
effect is related to the four-component spinor structure, i.e., to the existence of two 
non-zero spinor components x' an d f' of the wave function. This will be shown in the 
following. 

Let us consider the stationary solutions of the Dirac equation by writing the spinor 
\1/ in the form 

y = e -iE/nthc\ f (9) 

| We prefer to use this name because there is no orbital motion for a free particle 
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where E is the total (kinetic plus rest) energy of the fermion. Inserting this expression 
into (0) we get two equations for the spinors cp and x 

(E — mc 2 )x — o ■ r (f ■ p-\ — L • a)c<p (10) 

r 

i 



(E + mc 2 )<p — a • r (f ■ p-\ — L ■ a) cx- (11) 

r 

From these equations the L — S coupling of the spinors <p and x becomes apparent. 
Let us consider solutions with spherical symmetry of these equations. In Appendix A 
a derivation slightly different from the one used in most textbooks is presented. It is 
shown that the spinors can be written as products of a radial and an angular function 

as 

X = iG ji (r)$ jtm (9,<l>) , (12) 
(p = -F jei (r)<f> jem (e, (f)) = F je ,(r)a ■ f <S> j£m , (13) 

where £' = i — k/\k\, k being a non-zero quantum number which has a different sign 
according to the way the spin couples to the orbital angular momentum (see ( A.3 )). 



Since i' ^ i the whole spinor ^ is not an eigenstate of the orbital angular momentum 
operator L 2 . The good quantum numbers are j (total angular momentum quantum 
number), s — |, m (see Appendix A) and parity. This is due to the L — S term 
mentioned above. 

It is interesting to look at the non- relativistic limit of the equation (|ll|)- If we 
divide it by mc 2 , we obtain 

— + 1W> = — cr-px, 14 

,mr / mc 

using the fact that a ■ r(r ■ p+ l -L ■ a) c = a ■ pc. In the non-relativistic limit, the linear 
momenta of the dominant plane-wave components of x (obtained through a Fourier 
decomposition) are much smaller than mc, which implies, from (0), that if disappears 
in that limit. Since the angular part of tp contains only geometrical information, we 
can conclude (see equation (|T3"D) that Fjgi vanishes in the non-relativistic limit and one 
recovers the two-component spinor description of a spin-| particle. 

Interestingly enough, in the ultra-relativistic limit, where E + mc 2 ~ E — mc 2 ~ E, 
we can again recover the two-component description, since in this case we can choose 
the spinors x an d V 9 to be eigenstates of the helicity operator 

where p = p/\p\, with eigenvalues ±1, as can be seen from equations fllOD and ([TT|). This 
implies that x = an d therefore we may construct two two-component spinors for 
each value of the helicity (see, for instance, Itzykson and Zuber 0). 

If we differentiate once the coupled first order differential equations for Gj£ and Fj? 
derived in Appendix A (equations ( |A.11| ) and (|A.12| )), one gets 

d 2 G E + 2dG E _ £{£+!) Ct + E 2 -m 2 c* c ^ = Q 
dr 2 r dr r 2 5 (he) 2 ~* 
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d^F> + 2 dF> _ + + g 2 -mV fi = Q 

dr 2 r dr r 2 5 (/ic) 2 5 

Notice that, although in each equation only I or appear explicitly, the radial functions 
depend also on j through the energy E. These are the differential equations which have 
to be solved in order to get the radial functions Gji and F^. 

In the non-relativistic limit, since Fjg vanishes, i is again a good quantum number, 
i.e., the L — S coupling disappears. Moreover, since in this case Gj£ only depends on £, 
we can construct the standard non-relativistic solution taking the linear combination 

G e ^2(£ mi ; ]- m s \ j m)$ i£m = GeY eme Xm s , (17) 

jm 

where we dropped the index j in the radial function and used the definition (|A. 1|) of 
§jim and an orthogonality property of the Clebsch-Gordan coefficients. 

The differential equations ( |A.11|) and (|A.12|) can be extended to include interactions 
with spherical external potentials V(r) and m(r), which are respectively a time 
component of a four-vector (affecting the energy) and a Lorentz scalar (affecting the 
mass). This is done by the replacements E — > E — V(r) and m — > m(r). 

3. Solution of the Dirac equation in an infinite spherical well 

In order to show numerically the relativistic nature of the L — S coupling described 
in the preceding section, we are going now to compute the positive energy solutions 
of the Dirac equation for an infinite spherical well. As we will show, we can go, in a 
natural way, from a relativistic to a non-relativistic situation by changing the radius of 
the potential. The boundary conditions at the wall of the potential provide a discrete 
energy spectrum which allows a clear picture of the non-relativistic limit. 

To solve the Dirac equation in such a potential, one has to avoid any complications 
due to the negative energy states when trying to localize a spin-^ particle within a 
distance of the order of its Compton wavelength h/(mc) or less (this is the case for 
confined relativistic particles), one example of which is the Klein paradox (see, for 
instance, [|J). In other words, we want to retain the role of the Dirac equation as a one- 
particle equation in the presence of a infinite external potential. This is accomplished 
by defining a Lorentz scalar potential, i.e., a mass-like potential, having the form 

m r < R 

(18) 

oo r > R 

where m is the mass of the particle. The effect of this potential is to prevent the 
particle from propagating outside the well, meaning that its wave function is identically 
zero there. Inside, it behaves as a free particle of mass m. A potential like (|T8|), usually 
with m — 0, has been used to describe confined quarks as constituents of the nucleon 
in the MIT bag model (see, e.g., [[| for a review of this and related models). 

The boundary condition for the wave function at the boundary (r = R) cannot 
be obtained by requiring its continuity, since, being the Dirac equation a first-order 



m(r) 
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differential equation, the potential (JT3) implies that there is an infinite jump in the 
derivative of \I> (i.e., in the radial derivatives of G and F) when the boundary of the 
well is crossed. This jump obviously would not exist if \l/ were continuous. Another and 
most natural alternative is to demand that the probability current flux at the boundary 
is zero. As it is shown in ||, this is also a necessary condition to assure the hermiticity 
of the kinetic part of the Dirac hamiltonian within the well. This can be achieved by 
the condition 

-i(3a-fV = # at r = R. (19) 

In fact, if one multiplies this equation on the left by W(3 and its hermitian conjugate 
on the right by (3^! one gets -i^a ■ f\P = ^(3^/ and -i^a ■ = at r = R. 

These two equations imply that and \E^<3 • f\l/ are zero at r = R. 

The expression ■ f^> can also be written as j ■ r/c, where j = ^a^c is the 
probability current density for the particle described by the wave function Instead 
of the current flux we can look at the value of W (3^ at the boundary: indeed, since \l/ 
is zero for r > R, we may as well summarize the effect of the boundary condition (|19|) 
by saying that ^/3^f is continuous for any value of r. 

Having established the boundary condition, we proceed now to compute the radial 
functions. This is done in Appendix C. The full spinor \1/ reads 

where A is determined from normalization. In order to obtain the energy spectrum, we 
apply the boundary condition (|I9D to the spinor (p0|) . This gives rise to an equation 
relating the two radial functions (see Appendix C) 

j t (X )= J —MX 21 
I ft I V fcj + ror 

where X = VE 2 - m 2 c A /(hc) R. It can be written, in a more convenient way, in terms 
of the scaled quantities y = (E — mc 2 )/{mc 2 ) and xr = R/L , with L = h/(mc). 
These are the kinetic energy in units of mc 2 and the well radius in units of the Compton 
wavelength, respectively. We get then 

je(x R y/y 2 + 2y) = - * /IZI j e (x R ^y 2 + 2y) . (22) 
v |ft| y y + 2 v 

This equation is solved numerically for y as a function of xr for a given set of £, £' and 
ft. The results are presented in Figure 1. We plot the first values of y up to £ = 5 
for three values of xr. The energy levels are labeled in standard spectroscopic notation 
n£j, where n denotes the nth solution for a given set of £ and j. For xr = 100 the 
non-relativistic results, using the notation n£, are also presented. The non-relativistic 
spectrum is obtained by solving the Schrodinger equation for a particle of mass m in an 
infinite spherical potential well of radius R = lOOLo. The solutions can be found, for 
instance, in the quantum mechanics textbook of Landau jnj (in this case, there is no 
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Figure 1. The first 17 scaled kinetic energy levels obtained by solving equation 
for values of xr equal to 1, 10 and 100. In the last plot the first 10 scaled kinetic energy 
levels obtained by solving the Schrodinger equation for an infinite spherical potential 
well of radius R = IOOLq ar e also shown. 



spin-orbit coupling of the type (P because the potential is zero inside the well). The 
radial functions are spherical Bessel functions subject to the boundary condition 

MkR)=0, k = ^j^, (23) 

where is the kinetic energy of the particle. Notice that, in the non- relativistic limit, 
y <C 1, equation (|22"D reduces to ( p3|) since 

x R y/y 2 + 2y~ x R ^Ty = R , (24) 

and the factor \Jy/{y + 2) goes to zero in this limit. 

Analyzing Figure 1, we see that, as the radius of the well increases, the energy levels 
with the same I start grouping until they become degenerate and almost identical to 
the corresponding non-relativistic values. This effect is more pronounced for the states 
with higher t (notice the behaviour of the I/19/2 and l/tn/2 states). So we can conclude 
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that going from a radius R = L to a radius i? = 100L the L — S coupling effect fades 
away and j is no longer needed to classify the eigenstates of the system, and instead the 
orbital momentum quantum number £ emerges as the relevant quantum number. Since 
the boundary condition ( |22| ) effectively imposes a (maximum) value for the wavelength 
of the wave function and thereby a (minimum) value for the energy through the De 
Broglie relation, increasing the radius of the well amounts to decreasing the energy until 
we reach non-relativistic values for R = 100L . Notice that for the higher levels, for this 
value of R, even though there is not a perfect match with the non-relativistic energy 
values, the vanishing of the L — S coupling is a fact. The crucial scale here is the 
Compton wavelength Lq = h/(mc), determining the relativistic nature of the solution 
through the well radius. 

In summary, we have showed numerically the relativistic nature of the L—S coupling 
in the Dirac equation by computing its solutions for a particle with mass m in an infinite 
spherical potential well of radius R and making R sufficiently big as to produce non- 
relativistic solutions. 
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Appendix A. 

In this Appendix the radial equations for the Dirac equation are derived. We first write 



(f and x m equations fllPD and (|TTD as products of a radial function and a function 
of the angular coordinates 9 and <fi. To be able to get ordinary differential equations 
for the radial functions, the angular function must be an eigenstate of the operator 
L ■ a = (J 2 — L 2 — S 2 ) I '% (where J — L + S denotes the total angular momentum), which 
acts only on the angular coordinates. Accordingly, the angular function, $j£ m , reads 

$jtm(0,4>)= J2 J2 ^m t \ - m s \jm)Y £me (6,(j))xm s (A.I) 

mi_=-lm a =-\ 

where Ye m (8,<j)) is the spherical harmonic with quantum numbers i and m, Xm s the 
two-component spinors 



1 



and (Img; \ m s \jm) is a Clebsch-Gordan coefficient. The wave function ( |A.1|) is an 
eigenstate of J 2 , L 2 , S 2 and J z with eigenvalues h 2 j(j + 1), h 2 £(£ + 1), h 2 \(\ + 1) = h 2 \ 
and km respectively. Therefore we have 

L ■ 3 $ jim = h\j{j + 1) + £(£ + 1) - 3/4]% m = 

= -h(l + K)$ jtm (A.2) 
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with 



K 



(A.3) 



j = e + 

{ t= j + \ 3=1- _ 

For a fixed j, the quantum number k takes into account the two different possibilities 
for £, namely £ = j ± |, by just changing its sign. It also satisfies the equality 
k(k + 1) = £(£ + 1) for a certain I. Thus k can be considered as an alternative 
quantum number for the wave function Qje m replacing I. The corresponding operator 
is — (H + L ■ a). Note that for i = only one value of k is defined (-1). Wave functions 
with a fixed j but different £'s have opposite parity since £ = j ± ~ and parity is given 
by (—1)^). Using standard notation, and x are then written as 

X = iG jt {r)^ jem {e i <f>) , (A.4) 
(p = -Fj t ,(r)<$> jVm {6, <p) . (A.5) 

The quantum number £' of the lower component tp can be found by applying the operator 
a ■ r(r ■ p+ X -L ■ a)c to x ( see equation ([TT|)), giving 

i -> d i -> 

a ■ r(r ■ p+ -L ■ a)cx = ca ■ f(-ih— + -L ■ B)\G jt {r) <& jlm 



he 



dG 



, +(l + «)^ 
dr r 



<7 • f $«, 



jlm 



(A.6) 



The effect of a ■ r over $j^ m can be computed using the tensor properties of a and f (see 
Appendix B), yielding 



a ■ f $ j£m = 
where £' is given by 

+ 1 



jt'm j 



(A.7) 



_ 1 
2' 



if 3 
-1 if j 

Note that £' is related to £ by £' — £ — k/\k\. If we define the operator 

(-(h + L-a) \ 
^ h + L-a) ' 

^ will be an eigenstate of K with eigenvalue k. Thus k is also a good quantum number. 
From (|TT|) and ( |A.6|) we can write <p in (|A.5| ) in the form 

ip = F je ,(r) a ■r<b jlm {Q 



(A. 



(A.9) 



(A.10) 



The radial functions Gji(r) and Fjp(r) satisfy the coupled differential equations 
(see equations (|T0| ) and (JTTj) ) 



— mc 2 )Gj£ = —he 
(E + mc 2 )F jfJ = he 



dFj/r , .. Fjfi 
" +(1 + «') — 



dr 
dGj£ 



r 



, +(l + «)^ 

dr r 



(A.ll) 
(A.12) 



where is related to £' in the same way as in ( |A.3| ) (giving the relation k' = —k) and 
the relation a ■ r<&je' m = —<&je m was use d (note that (a ■ f) 2 = I). 
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Appendix B. 

In this Appendix we will derive expression ( |A.7| ) by calculation the matrix element 

®ri>m>°-rQjtm, (B.l) 

where '+' stands for hermitian conjugate. Since both a and f are vector operators 
(irreducible tensor operators of rank 1) use can use a general theorem for the matrix 
element of a scalar product of commuting tensor operators bet ween eigenstates of 
angular momentum. Using the notation and conventions of Edmonds |7| we have 

^a-fV^-lfH/^'lj \ i} (lll^lll) (?¥¥) ■ (B.2) 

Using the conventions of Edmonds, the reduced matrix elements can be evaluated, such 
that O is 

i / f i 1 ^ ^ 

*T 7f . £ if\ . . —I 1 ^+o+.?x X atno i 1W1 n. inl tf'n\ J J 



<5>} e , m ,a-r^ em = (-l) £+ 2+^ y 5 wv /6(2£ + l)(10;£0|£ , 0)| J 1 J ^ ( B . 3 ) 

where we used the fact that (10; £ | £' 0) is non-zero only for = £± 1. The 6-j symbol 
is different from zero only for j = ± |. Since we have also j = £ ± ~, we have two 
possibilities for a fixed £: 

1) i = / + i =► * = e + i ; j = '-^ (b.4) 

2) j=l-\ =► * = <-l ; J = ^' + ^ (B.5) 

Inserting the values of the Clebsch-Gordan coefficient (10; £ | £' 0) and of the 6-j 
symbol into (|B.3|) we get, for both cases, 



4«'m' ^ " ^ ^jtm = Sjj'Srnm' , (B-6) 

where £ and are related by ( p.4| ) and ( |B.5| ). Since the spinors $je m form a complete 
orthonormal set this equation implies ( |A.7|) . 



Appendix C. 

In this Appendix we obtain the spinor which is the solution of the Dirac equation with 
the infinite spherical potential (0). To compute the radial functions Gje and Fji> inside 



the well, we first look at equations (15) and (|T^) and make the change of variable 



x = E 2 — m 2 c 4 / (he) r. In this way, we get equations of the form 
d 2 /i 2d/, / + 



where / and // stand for £, £' and Gji, Fj#, respectively. The solutions of equation ( |U.1| ) 
which are regular at the origin are the spherical Bessel functions of the first kind, ji(x) 
(see, for instance, Abramowitz and Stegun ||). Since these solutions are determined up 
to an arbitrary multiplicative constant, the radial functions are 

G jt = Ajt(x) (C.2) 
F 3e = B 3e (x) , (C.3) 



On the relativistic L — S coupling 1 1 

where A are B are constants. We can use one of the equations (|A.ll ) or ( |A.12| ) and the 
recurrence relations of the functions ji(x) (see H]) to find the following relation: 



B = (C.4) 
\k\ V -c/ + mc z 

The complete spinor \l/ then reads 

• ■ (r 8 4 t) - Ae-"» ( 'M^^r^jOA) \ 
v 3Km \r,v,<p,t) - Ae | / g - me 2 . ^ //i jl\ • 

ft: — o lt>\ =- r sP^/r 



Applying the boundary condition (|T^) to this spinor leads to 



where X = \/E 2 — m 2 c 4 / (he) R and the relation ( |A.7| ) and its inverse were used. This 
equality implies 



MX) = ~J^rr^MX) (C.7) 
\k\\L + mc z 
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